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ABSTRACT 


Like most design problems , kinematic synthesis of planar linkages involves a number of 
trials to reach an acceptable or an optimal solution . Both graphical and analytical metho- 
ds are available for the synthesis of four-link planar mechanisms . Independent of the me- 
thod used , a solution arrived at always needs to be verified through animation . Consequ- 
ently , a versatile software incorporating graphical and analytical methods and capable of 
analysis , synthesis and animation will be very useful to a mechanism designer . 

This thesis , as a part of a D. S. T. ( Department of Science and Technology ) 
sponsored project , attempts to develop some specific problems , the solution of which 
require special graphical techniques or overall optimization problem . The Borland C H ' 
language with graphic facility has been used because of its certain special features , which 
are eminently suitable for the problems under consideration . The software informs the 
user the kind of problems it can handle , asks for the input data ( including designer ’ s 
choice ) , verifies the feasibility of the existence of solution , obtains the optimal solution 
and demonstrates the performance of the obtained solution . 
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1. INTRODUCTION 


1.1 INTRODUCTION 

Planar linkages are extensively used in various machineries in the field of agriculture, 
printing, textile and many others. A successful design of such a machinery very often 
requires, as a first step, an optimal kinematic design of various planar linkages used in it. 

The kinematic design of a linkage starts with the number synthesis, i.e., with the 
decision on the number of links and kinematic pairs to be used for attaining the desired 
objective. However, a well-developed theory exists only for designing four-link 
mechanisms. The theory also points out, in some cases, the limitations of such four-link 
mechanisms. Only when the job in hand can not be achieved by a four-link mechanism, 
one needs to go for more complicated ( e.g., six-link or eight-link ) mechanisms. Even for 
designing such higher order linkages, very often the theory of four-link mechanisms 
serves as the starting point. 

Over the last hundred years, various graphical and analytical techniques have been 
developed for the dimensional synthesis ( i.e., the determination of the actual kinematic 
dimensions ) of four-link planar mechanisms. The suitable choice of a method depends on 
the type of problem. In a four-link mechanism, the links connected to the fixed link ( 
ground or frame of reference ) are referred to as the input and output links, whereas the 
link not directly connected to the fixed link is called the floating link ( e.g., the coupler, 
the connecting rod ). 

Depending on the required kinematic characteristics to be satisfied by the 
designed linkage, dimensional synthesis problems can be broadly classified as 
[i] Motion Generation - It is defined as the control of a line segment in a plane ( on the 
floating link ) such that it assumes some prescribed set of sequential positions. The 
guidance may or may not be coordinated with the input motion. 
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[ii] Path Generation - If a tracer point on the floating link has to traverse through some 
specified points ( on a specified path ), then such a problem is classified as a problem of 
path generation. The generation of a prescribed path may or may not be coordinated with 
the input motion. No attempt is made to control the orientation of the link which contains 
the tracer point. 

[iii] Function Generation - In this problem, the movements of the output and input links 
are coordinated according to a prescribed function [ say, y = f (x) ]. The change in the 
input movement corresponds to the change in the independent variable and that in the 
output movement to the change in the dependent variable. A function generator is 
conceptually a “ black box ” which delivers some predictable output in response to a 
known input. 

[iv] Dead-Center Problems - These types of problems generally occur where the quick- 
return action of the output link for a constant input speed becomes important. Here; the 
linkage dimensions are determined for fulfilling the prescribed dead-center configura- 
tions. 

In a real-life problem of path generation or function generation, it is generally not 
possible to design a mechanism which achieve the specified characteristics over the 
desired range of movements. One can only attempt for an approximate solution. There are 
two approaches to achieve an acceptable level of approximation. 

[a] Accuracy Point Approach - In this approach, the generated characteristics of the 
designed mechanism exactly match with the desired characteristics only at some discrete 
positions within the specified range of movement. But at remaining configurations there 
exists some difference between the desired and generated characteristics which is called 
the structural error. One should select those precision points within the specified range 
of movement for which the maximum structural error becomes as small as possible. As a 
first approximation, Chebyshev’s accuracy points [l) are generally chosen. With this 
approach, we can get a better solution if we take more number of precision points within 
the specified range of movement, but the number of precision points is limited by the 
number of available design parameters. 
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fb] Optimization Approach - In contrast to the above approach, here, the structural 
error, which may not be zero anywhere, is minimized in an overall sense in the entire 
range. It should be noted that, one can choose as many design positions as he / she wants. 
The number of design positions is independent of the number of available design 
parameters. 

The accuracy point approach can be adopted while applying either a graphical or 
an analytical method, but the optimization approach is amenable to only analytical 
methods. 

Finally, independent of the approach or technique adopted, a designed linkage has 
always to be checked for branch and'order defects by checking through the entire range of 
movement. Moreover, there are always some assumed parameters and an optimal solution 
can be attained only after a number of trials. 

All the points mentioned above make it amply clear that the dimensional synthesis 
of four-link mechanisms is an ideal problem to be tackled by a computer software with 
animation and graphic display. Commercial softwares are, of course, available, which are 
mostly suitable for kinematic and dynamic analysis. In this thesis, an attempt has been 
made, as a part of a D.S.T. ( Department of Science and Technology ) sponsored project, 
to develop a software which encompasses certain special graphical methods and 
optimisation approaches to solve a variety of dimensional synthesis problems of four-link 
planar mechanisms. 

1 . 2 OBJECTIVE AND SCOPE OF THE PRESENT WORK 

The objective of this work is to develop a computer software which can synthesize and 
animate four-link mechanisms for different types of problems, which are 

• Dead-center problems for crank-slider and crank-rocker mechanisms : Here 

[i] the quick-return ratio, 

[ii] the stroke-length ( in the case of 3R-1P ) or the rocker swing angle ( in the 

case of crank-rocker ) and 

[iii] one of the link lengths or offset 
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have to be given as inputs. This kind of problem has to be solved while 
designing quick-return mechanisms. 

• Generation of coupler curves : To draw a segmented coupler curve corres- 
ponding to a specified coupler point on the floating link. Each segment 
corresponds to a specified crank-rotation. Thus the velocity of the coupler 
point at different positions can be judged by observing the length of each 
segment. Larger the length, greater is the velocity of the coupler point at that 
position. 

• Path generation by optimization approach : In this case, the coordinates of 
several points on the specified path, coordinated with input angles will be 
required as inputs. In some cases, where the shape of the path along with the 
structural errors at some specified points are important, this approach is 
usually necessary. 

• Function generation by optimization approach : Here, either a function or a 
set of coordinated input-output angles are prescribed. This is generally used in 
the synthesis of function generators. 

• Two cusps problems : Here, a crank-rocker has to be designed so that a 
coupler point curve has two cusps at prescribed locations. Further, the crank 
rotation as the coupler point moves from one cusp to another is also prescrib- 
ed. This type of problem arises in several processes ( e.g. in the stamping and 
feeding processes, where a coupler point has to reach at its destination and 
return back without impact and shock, noiseless rocker-stop used in type- 
writers, recording machine in which a recording pencil writes only between 
two limiting points etc. ). 

• To find the coupler curve cognates : In the case of a path generation problem, 
a designed mechanism may not be accepted if its fixed pivots are located 
outside the available space. Then, in such a situation, its cognate generating 
the same coupler curve may be accepted. Similarly, if the initially designed 
mechanism turns out to be a Grashofian double-rocker, then its crank-rocker 
cognate may be an acceptable design. 
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• Problems with specified kinematic variables ( angular velocities and 
accelerations ) : To design a 4R linkage so that at a particular configuration, 
the prescribed values of 

[i] angular displacements, velocities and accelerations of the input and output 
links or 

[ii] angular velocities and accelerations of the input, floating and output links are 
achieved. 

To meet the computational and animation requirements for tackling the 
above classes of problems, we have selected C ++ ( Borland, with graphics environment ) 
as the programming language because of its special characteristics ( e.g. classes, objects, 
inheritance, overloading operator, data hiding etc. ). The classes and some special 
member functions used in the present work are listed in the appendix. 



2. DEAD - CENTER PROBLEMS 


2. 1 SLIDER - CRANK MECHANISM 

Let us first consider the problem of synthesizing a slider-crank mechanism with 
prescribed stroke-length and quick-return ratio. 

Referring to Fig 2. 1 , suppose 

B[ and Bo represent the extreme positions of the slider 4 ( when the crank 2 
and the connecting rod 3 become collinear ), 

H = stroke-length of the slider = B1B2, 

Qrr = quick-retum ratio, 

0 2* = An angle through which the crank rotates when 
the slider goes from B] to B 2 
or 6 2 * - 2 tz { ( Qrr ) / ( 1 + Qrr) }, 

02 **= ZB, 0 2 B 2 = ( 0 2 * -it ) 
and - ZO2B1B2. 

If only H and Qrr ( or 0 2** ) are given, then we can have infinite solutions. It can 
be understood easily by using the property of a circle which states that “ a chord of a 
circle subtends equal angle at all the circumferencial points which are on the same side of 
the chord So O2 can be placed anywhere on the circumference of a circle ( CR2 ) 
passing through Bi, B2 and O2. Thus, for a unique design, besides H and 02* ( i.e., the 
quick-return ratio ), one more kinematic dimension, viz., L2 ( crank length ) or L3 ( 
connecting rod length ) or E ( offset ), must be either prescribed or assumed. Depending 
on which one of the three dimensions is given ( or assumed ), three methods of 
determining the other two are discussed below. 

Case 1 . H, 0 2 * and E are prescribed, L 2 and L 3 are to be determined : 


From A ’s 0 2 NBi and 0 2 NB 2 , 
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E = (L 3 -L 2 )sin( p + 0 2 ** ) = ( L 3 + L 2 ) sin ( p ) [2.1.1 ] 

and 

H = (L 3 + L 2 )cos( p ) - ( L 3 - L 2 ) cos ( p + e 2 **). [2.1.2] 

From Equations 2 . 1.1 and 2. 1 .2, 

sin ( p + 0 2 **) sin ( p ) = ( E/H) sin ( P + 0 2 **- p ) 
or 

cos (2 p + 02 **) = cos ( 0 2 ** ) - ( 2 E / H ) sin ( p + 0 2 **- p ) 
or 


Now 


and 


P = [ cos ’ 1 { cos ( 0 2 ** ) - ( 2 E / H ) sin ( 0 2 **) } - 0 2 ** ] / 2. 

[2.1.3] 

0 2 Bi = ( L 3 + L 2 ) = E / sin ( p ) [2.1.4] 


0 2 B 2 = ( L 3 - L 2 ) = E / sin ( p + 0 2 **). [2.1.5] 

Therefore, 

L 3 = (0 2 B, +0 2 B 2 )/2 and L 2 = ( 0 2 B, - 0 2 B 2 ) /2. [2.1.6] 

Case 2. H, 0 2 * and L 3 are prescribed, L 2 and E are to be determined : 


From A 0 2 B 2 Bi, 

cos ( 0 2 ** ) = { ( L 3 + L 2 ) 2 +( L 3 - L 2 ) 2 - H" } / { 2 ( L 3 + L 2 ) ( L 3 - L 2 ) } 
or 

cos( 02 **)={L 3 2 + L 2 2 -(H 2 /2) }/(L 3 2 -L 2 2 ), [2.1.7] 

which gives 

L 2 = /"[ [ ( H 2 / 2 ) - L 3 2 { 1 - cos ( 0 2 ** ) } ] / { 1+ cos (0 2 **)} ]• 


[ 2 . 1 . 8 ] 
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From Equation 2. 1 . 1 , we get 

( L 3 - L 2 ) { sin ( j3 ) cos ( 6 2 ** ) + cos ( j3 ) sin ( 6 2 ** ) } 

= ( U + L 2 ) sin ( (5 ) 
or 

tan ( (3 ) = { ( L 3 - L 2 ) sin ( 0 2 ** ) } / { ( L 3 + L 2 ) - ( L 3 - L 2 ) cos ( 0 2 ** ) } 
or 

P = tan 1 [ { ( L 3 -L 2 ) sin ( 0 2 ** ) } / { ( L 3 + L 2 ) - (L 3 - L 2 ) cos ( 0 2 ** ) } ]. 

[2.1.9] 

Hence 

E = (L 3 + L 2 )sin( p ). [2.1.10] 

Case 3. H, 6 2 * and L/> are prescribed, L 3 and E are to be determined : 

From Equation 2. 1 .7, 

L 3 = a T[ [ ( H 2 / 2 ) - L 2 2 { 1 + cos ( 0 o** ) } ] / { 1 - cos ( 0 2 ** ) } ]• 

[2.1.11] 

Then using Equations 2.1.9 and 2. 1 .10 we can get the offset E. 

Optimal solution : As we have seen above, that besides H and Qrr (or 6 2 * ), One has to 
give either L 2 or L 3 or E to obtain a unique solution for a 3R-1P mechanism. If only H 
and Qrr ( or 6> 2 * ) are available then one can also design a unique slider-crank mechani- 
sm for which the minimum transmission angle is maximum. 

For smooth functioning of a mechanism in real life, we should ensure that, besides 
satisfying the kinematic requirements, the mechanism moves freely. Neglecting the 
effects of the gravity, inertia and friction all binary links become two force members. 
Consequently, the force in any such link is along the link vector and the force analysis 
becomes simple. To judge the free running quality of a mechanism, the transmission 
angle is considered as one of the most important parameter. 
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For a slider-crank mechanism the transmission angle is defined as the absolute 
acute angle between the connecting rod and the normal to the slider movement ( Fig - 2.2 
). The minimum transmission angle [1] for this mechanism is given by 

Atmin = cos- 1 [(L 2 + E)/L 3 ]. [2.1.12] 


Referring to Fig 2.1, 

from A O2B1B2, we observe that 

L 3 -L 2 = { H /sin ( 0 2 **) } sin (p) 

and 

L 3 + L 2 = { H / sin ( 6 2 ** ) } sin ( + 6 2 ** ). 

Therefore, 

l 2=[ 2 sin 09 2 ]{ sin (^ + 0 2 **)-sin(/ 3 )} [2.1.13] 

and 

U=[ ](sin(/J + e 2 **) + sin(^)J. [2.1.14] 

Again, from A 0 2 B[N, we get 
E = ( L 3 + L 2 ) sin ( p ) 
or 

E = [ , — ] sin ( P + 0 2 **) sin ( p). [2.1.15] 

L sm(02**) 

To maximize the minimum transmission angle jl min ( given by Equation 2.1.12 ), 
one must minimize the value of [ ( L2 + E ) / L 3 ] which is a function of the variable fi. 
Using Equations 2.1.13, 2.1.14 and 2. 1.15, we get 

L 2 + E r 2 si n ( f3 ){ 1-sin ( + %**)} -i roii -51 

L 3 " sin ( /? + 0 2 * * ) + sin ( ) J ' 1 
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From Fig 2.1, it is obvious that 0 2 can not be above the line passing through Bi 
and Bo ; and 0 2 Bi must not be less than 0 2 B 2 ( since B| and B 2 are the outer and inner 
dead-centers, respectively ). This implies that ft varies from 0 to (k - 0 2 ** ) / 2. 

Thus, we have a function [ ( L 2 + E ) / L 3 ] of a single variable /?. The variable /? 
varies from 0 to ( n - 0 2 ** ) / 2 . As we know that, to minimize a single variable func- 
tion for a given range of the variable, one can use the golden section search method [2] . 
For the present problem, using this optimization method, we get /? opt corresponding to the 
minimum value of [ ( L 2 + E ) / L3 ] and hence L 2 , L3, E and fi m j n ( from Equations 
2.1.13, 2.1.14, 2.1.15 and 2.1.12 respectively ) are known. 

Flow chart with illustrations - Using the computational flow chart ( Chart 2. 1 ), one can 
easily understand the steps used in the program for dead-center problems of a slider-crank 
mechanism. It also helps the designer to run the program. After computing the design 
parameters of the mechanism, the software itself checks the existence of the mechanism. 
If the mechanism exists, then DRAW_3R1P( ) can be used to animate the slider-crank 
mechanism. For verification of the results, one should opt for dead-center configurations 
which display all the dead-center configurations along with the stroke-length. 

To illustrate the synthesis of a slider-crank mechanism ( case 1 ), we have 
considered a problem with H = 6.0 unit, Qrr =' 1 .4 and E = 1 .0 unit. 

For this problem, the final design yields 
L 2 = 2.8629 unit, 

L 3 = 4.489 unit 
and fi min = 30.73°. 

To verify the result, one can see its dead-center configurations shown in Plot 2.1. 

For an optimal design problem, taking H = 6.0 unit and Qrr = 1.4, we design a 
slider-crank mechanism for which ' 

L 2 = 2.7307 unit, 

L 3 = 5.5222 unit. 



E= 1.9198 unit 
and /! min = 32.63°. 

Its dead-center configurations are shown in Plot 2 . 2 . 

2. 2 CRANK - ROCKER MECHANISM 

Now let us consider the synthesis of a crank-rocker mechanism for dead-center problems. 
From Fig 2.3, O 2 A 1 B 1 O 4 ( shown in Fig 2.3 ) is a crank-rocker mechanism. The rocker 
O4B1 swings by an angle 04 * when the crank CFA] rotates by 02 * ( CCW ) and 
correspondingly the moving hinge of the follower traverses from Bi ( outer-dead center ) 
to Bt ( inner-dead center ). The link lengths ( L|, L 2 , L 3 and L 4 ) are indicated in the 
figure. 

Here B 1 B 2 = 2 L 4 sin ( 0 4 * / 2 ) and the offset E can be taken as the perpendicular 
distance of line B 1 B 2 from 02 - So if the length B 1 B 2 ( = the stroke-length ), Qrr (or 0o* 
) and one of the lengths ( Li or L 2 or L 3 or E ) are given, then one can synthesize a crank- 
rocker mechanism in the same way as was done for the slider-crank mechanism. 

But in real life, generally Qrr ( or 0 2 * ), the rocker swing angle 0 4 * and the fixed 
length Li are given to synthesize a crank-rocker mechanism. In such a situation, the 
solution is obtained by the Alt’s construction m . The Alt’s construction gives infinite 
solutions, but one must' select the mechanism for which the minimum transmission angle 
is maximum considering the free running quality of the mechanism. 

The transmission angle of a 4R mechanism is defined as the absolute value of the 
acute angle between the output link and the coupler ( Fig 2.4 ). Suppose a torque T 2 is 
applied at link 2, then the generated torque T 4 (on link 4 ) will be directly proportional to 
the force F t [ = F sin( ) ] where /I is the transmission angle and F is the produced axial 

force in the coupler. The radial component F r provides only tension or compression in 
link 4. For good performance of a mechanism it is suggested that ju min > 30° - 35 °. 

The transmission angle [1] /i (for the crank-rocker ) is given by 


/.l = cos " 1 [ { ( L 3 2 + L 4 2 ) - ( Lf + L 2 - ) - 2 Li L 2 cos 02 } / ( 2 L 3 L 4 ) ]. [ 2.2.1 ] 
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If 6 2 * > n [ or ( L 3 2 + L 4 2 ) > ( L, 2 + L 2 2 ) ] then 


d- min — cos [ { ( L 3 ~ + L 4 ~ ) - ( Lj“ + L 2 " ) + 2 L[ L 2 } / ( 2 L 3 L 4 ) ] 


and if 0 2 * < k [ or ( L 3 2 + L 4 2 ) < ( L, 2 + L 2 2 ) ] then 


[ 2 . 2.2 ] 


A I min = K - COS' 1 [ { ( L-f + L 4 2 ) - ( Lf + L 2 2 ) - 2 L, L 2 } / ( 2 L 3 L 4 ) ]. 

[ 2.2.3 ] 

To get the optimum solution, the designer generally uses Volmer’s nomogram ( 
Plot 2.4 ). To design a crank-rocker mechanism for prescribed dead-center configurations 
on a computer, one must have to feed the data from Volmer’s nomogram. But here, we 
have derived a single variable function ( say ji m - m ) which is again optimized by the 
golden section search method to obtain the corresponding unknown design parameters ( 
L 2 , L 3 and L 3 ). Finally we verified the results with the Volmer’s nomogram. 

Using the Alt’s constructions, the unknown design parameters of a crank-rocker 
mechanism for various situations are presented below. In all the cases 0 2 * and 0 4 * are 
measured CCW from the outer to the inner dead-center configurations. 

Let 

e 2 **=e 2 *- 7 t and 624 = 02 - 64 . [2.2.4] 

Case 1 : [ 024* < 7T and 0 2 * ^ n ] 

Referring to Fig 2.5, 0 2 and 0 4 are the fixed pivots such that 0 2 0 4 = Lj. At 0 2 , a 
line 0 2 E 2 is drawn at an angle of ( - 0 2 * / 2 ) to 0 2 0 4 , and another line 0 4 E 4 is drawn 
through 0 4 at an angle of ( - 0 4 * / 2 ) to 0 2 0 4 . These two lines intersects at F and the 
circle CR| is drawn with 0 2 F as the diameter. The mid-normal of 0 2 F intersects 0 4 F at 
C 2 . Another circle CR 2 , is drawn with its center at C 2 and radius as C 2 F. 
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In this case, since C2F is greater than ( if 9 2 * > n ) or equal to ( if 62* = n ) 
C2O4. The circle CR2 intersects O2O4 at S which coincides with O4 ( if 9 2* = it ) or lies 
outside ( if 62* > it ) O2O4 ( as shown in Fig 2.5 ). 

A line O4N is drawn through O4 at an angle of ( - 64* ) to O2O4 intersecting the 
circle CR2 at the point N. The point Bi can now be taken anywhere on the arc NS. The 
line O2B1 intersects CRi at A t . Thus, O2A1B1O4 is the required 4 R crank-rocker linkage 
at its outer dead-center configuration. All the unknowns and some of the assumed 
symbols are shown in the figure. 

Now to determine the unknown design parameters, referring to Fig 2 . 5 , 

let 

Ci = (x c) , y c i ) , ri = the radius of circle CRi, 

C2 = (x C 2, Yc 2 ) and x-i - the radius of circle CR2. 

From A O2FO4 


r, = 0 2 F/2 


r, = { Li sin ( 0 4 */2) } / { 2 sin ( 0 2 4 */ 2 ) } 

[2.2.5] 

[v Z 0 2 F0 4 = 024*/ 2], 


Xd =-n cos( 9 2* 12) 

[ 2.2.6 ] 

y ci =r ( sin ( 9 2* 12). 

[ 2.2.7 ] 


From A C 1 FC 2 

r 2 = FC 2 = ri sec ( 9 24 */ 2 ), 
from A O 2 FO 4 

X c2 = 2 X C | + T 2 cos ( 9 4* / 2 ) 

[v Z0 2 0 4 F = ( 9 4 *)/2 ] 

and 

y c2 = 2y c i -r 2 sin( 9 4 * 12 ). 


[ 2 . 2.8 ] 

[ 2.2.9 ] 

[ 2 . 2 . 10 ] 
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Let 

N = ( x n , y n ) and 0 4 N = l n , 

where 

x„ = Li -l n cos ( 0 4 *) 
and y n = l n sin ( 64 * ). 

Since N is on the circle CR2, we can write 
( X n - x c2 ) 2 + ( y n - y C 2 ) 2 = r 2 2 . 

Putting the values of x n and y n in the above equation 

l n 2 + { 2 ( x c2 - Li ) cos ( 9 4 * ) - 2 y c2 sin ( 6 4 * ) } l n + ( Lj 2 - 2 Li x c2 ) = 0. 

[2.2.13 ] 

Equation 2.2.13 gives two values of l n . Here, the positive value of l n is chosen as y n is 
positive. Then, from Equations 2.2.1 1 and 2.2.12, we can have x n and y n . 

Now 

tan fi „ = y n / x n 

or p n = tan _1 (y„/x„) . [2.2.14] 

For the point S, p s = 0. 

Thus, P can vary from P s to p n . 

Now let us consider a general point Bi on the circle CR 2 at an angle P ( which 
may vary from P s to P „ ). 

Suppose O 2 B 1 = L = L 2 +L 3 . 

Therefore, 

( L cos p - x C 2 ) 2 + ( L sin P - y C 2 ) 2 = r 2 2 
or 

L = 2 ( x c2 cos p + y C 2 sin /3 ) [ v x c2 2 + y c2 2 = r 2 2 ]. [ 2.2.15 ] 

Similarly the point Ai lies on the circle CRi at the same angle p , therefore. 


[ 2 . 2.11 ] 
[ 2 . 2 . 12 ] 


and on simplification, we get 


( L 2 cos p - Xd ) 2 + ( L 2 sin j9 - y c i ) 2 = rp 
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or 

L 2 = 2 ( x cl cos (3 + y c i sin /? ) [ v x c i 2 + y c i 2 = rj 2 ]. [2.2.16] 

Now 

L 3 = L-L 2 [2.2.17] 

and from A 0 2 Bi 04 

L4 = (L 2 + Li 2 -2L,Lcos p ) m . [2.2.18] 

Thus, for any angle ft (varying from j3 s to /3 n ) we can have the link lengths ( Li, 
L 2 , L 3 and L4 ) to design the crank-rocker mechanism. But we have already discussed that 
the mechanism for which fi m i n is maximum will be the best regarding the transmi-ssion 
quality. Now we have to maximize the minimum transmission angle ( given by Equation 
2 . 2.2 ) 

Mmin = cos-‘ [{ (L 3 2 + L 4 2 )-(L 1 2 + L 2 2 ) + 2L 1 L 2 ] / (2L3L4)]. 

[2.2.19] 

Since L 2 , L 3 and L4 are functions of /3 , hence jl mj n is also a function of the same 
variable /? (varying from [3 S to /3 n > To maximize we have used the golden 

section search method. As we know that this optimization method minimizes a single 
variable function in a given range of its variable. Therefore, in the present problem, we 
have chosen ( - fi m in ) ( since jU min is to be maximized ) as a function and j5 as a 

variable. 

After optimization, we can get j3 op t ( corresponding to ( jl min ) max ) and hence the 
desired unknown link lengths ( L 2 , L 3 and L4 ). 

Case 2 : [ 6 u* < k and 9 2 * < tc ] 

Referring to Fig 2.6, in this case, the point 0 4 lies outside the circle CR 2 because 
FC 2 < C 2 0 4 . Hence, the line 0 4 N (.at an angle - 0 4 * ) cuts the circle CR 2 at two points 
both with positive y coordinates. 

From Equation 2.2.13, one value of l n corresponds to the point N and the other to 
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the point S. Therefore, using Equations 2.2.1 1, 2.2.12 and 2.2.14 we can have and fi n . 
Further, from Equation 2.2.3, 

JU min = n - cos' 1 [ { ( L-f + L 4 2 ) - ( L, 2 + L 2 2 ) - 2 L, L 2 } / ( 2 U U ) ]. 

[ 2 . 2.20 ] 

Other Equations remain the same as in case 1. Again using the same optimization 
method, fi m \ n is maximized to get /3 op t ( corresponding to ( m i n ) max ) and hence the 
desired unknown link lengths ( La, L3 and L 4 ). 

Case 3 : [ 6 24* = n ] 

In this special situation ( Fig 2.7 ), the point C2 goes to infinity ( since Z 0 2 F0 4 = 
jt / 2 ) and the circle CR 2 degenerates into the straight line OaF. 

It is obvious that, OaF = FN since Z Oa0 4 F = Z F0 4 N = 64 * / 2. The point Bi 
can now be located anywhere on the line 0 2 F above N and the crank length ( = 0 2 F ) 
remains the same for all the solutions. 

From Fig 2.7, we observe that 

ft = k -( 62 * / 2), which is constant, [ 2.2.21 ] 

r, =L, sin ( 0 4 */2)/2 [2.2 22] 

and L 2 = 2 ri which is again constant. [ 2.2.23 ] 

Though the point S, in this case, can be at infinity, to achieve a compact linkage 
0 2 S is restricted to 4 0 2 F. Thus the variable L ( = 0 2 Bi ) is restricted in the range 2L2 ( = 
0 2 N ) < 0 2 B] < 4L 2 ( =0 2 S ). After optimization it has been found that it gives the same 
result as the Volmer’s nomogram. 


Now L3 = L - L 2 


[ 2.2.24 ] 
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and from A 0 4 BjD 

L 4 = { ( L sin /? ) 2 + ( L cos /? -L, ) 2 } 1/2 . [2.2.25] 

In this case, 

f* min = C0S‘‘[ { ( L 3 2 + L 4 2 ) - ( Li 2 + L 2 2 ) + 2 Li L 2 } / (2 L 3 L 4 ) ]. 

[ 2.2.26 ]. 

The golden section search method is again used to minimize ( - fl m j n ), but here the 

variable is L which varies from 2L 2 ( = L N ) to 4L 2 ( = L s ). Thus, we get L opt and hence 
L 3 and L 4 . 

Case 4 : [ 6 24 * > it ] 

This case is also similar to case 1 except for the change of location of the center 
C 2 of the circle CR 2 and as in case 2, the point S lies above the line 0 2 0 4 ( Fig 2.8 ). 
Therefore, 


r 2 = FC 2 = - n sec ( 6 24 * / 2 ) 

[ 2.2.27 ] 

and from A 0 2 F0 4 


x c2 = 2 x c i - r 2 cos ( 9 4 * / 2 ) 

[ 2.2.28 ] 

and 


y C 2 = 2 y c i + r 2 sin ( 0 4 */ 2). 

[ 2.2.29 ] 


Again Equation 2.2.13 gives two values of l n , one for the point N and the other for the 
point S. Then, from Equations 2.2.1 1, 2.2.12 and 2.2.14 we can have /3 S and /3„. 

Further, from Equation 2.2.2, 

fl m in = cos' 1 [ { ( L 3 2 + L 4 2 ) - ( Li 2 + L 2 2 ) + 2 Li L 2 } / (2L 3 L 4 )]. 

[ 2.2.30 ]. . 

Other Equations remain the same as in case 1. Then using the same optimization method 
fl min is maximized to get (3 op t [ corresponding to ( fl m m W ] and hence the desired 
unknown link lengths ( L 2 , L 3 and L 4 ). 

The following points regarding the Alt’s construction m may be noted. 
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[ i ] If a point Bi is taken outside the region NS, then the rocker has to cross the frame- 
line during its movement, which is not possible in crank-rocker mechanism. 

[ ii ] For a given value of 0 4 *, the minimum and maximum values of 0 2 * ( for a crank- 
rocker linkage to exist ) are given, respectively, by 

( 0 2* )min = ( K + 6 4 * ) / 2 and ( 6 2* )max = ( 3 n + 6 4 * ) / 2. 

[2.2.31 ] 

In these limiting conditions, the line 0 4 N in Alt’s construction becomes tangent to CR^. 

Flow chart with an illustration - The steps used in the programming are shown in 
Chart 2.2. The restrictions ( given by Equation 2.2.31 ) has been used to control the input 
data for achieving a feasible solution. Here, the role of the function GOLD ( type, p n , P s, 
xci, yci, Xc 2 j yc 2 ), used in the golden section search method is very important. The first 
variable of the function selects the type ( here, type = 1 for cases 1, 2, 4 and type = 2 for 
case 3 ) where as ( P n , P $ ) define the range of the variable and the last four variables 

determine the transmission angle through another function MEU (type, Xci, yci. Xc 2 ,.yc 2 > 
P ). In function GOLD ( ), a number r ( tau = 0.618) has been used, this is known as the 

golden number [2] . An accuracy controlling parameter e has also been considered in the 
same function, which controls the accuracy of (/i m in)max and p opt . We can get very 
accurate values just by controlling e , which is really an advantage over the conventional 
nomogram method. 

After knowing the link lengths, the function DRAW_4R( ) animates the 
mechanism. For verification of the results, one can use the option dead-center 
configuration. 

Now as an illustration let us consider a dead-center problem for which 
Li = 6.0 unit, 6 4 * = 60° and 0 2 * = 210°. 


The obtained final results are 
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( ft min ) 


max 


= 32.63°, 


P 0 p t = 46.57°, 

L 2 = 2.73 unit, 

L 3 = 5.519 unit 
and L 4 = 6.0 unit 

The dead-center configurations of this mechanism are shown in Plot 2.3. One can 
also verify ( /I m i n )max and jS op t from the Volmer’s nomogram ( shown in Plot 2.4 ). 




[Fig 2.1] 


O 2 A = L; 
AB = L 3 
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[ DEAD-CENTER PROBLEMS FOR 3R-1P MECHANISMS ] 



* Cdl -> [ I cos( 0 2 ** )-(2E/H)sin( 0 2 **) I < 1 ] 

* Cd2 -» [ ( H 2 / 2 ) > L 3 2 ( 1 - cos( 0 2 **))1 

* Cd3 [ ( H 2 / 2 ) > L 2 2 ( 1 +cos( 0 2 **))] 


[Chart 2.1 ] 
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[ DEAD-CENTER PROBLEMS FOR CRANK-ROCKER MECHANISMS ] 
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[ GOLDEN SECTION SEARCH METHOD ] 
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[Plot 2.1 ] 
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[ Plot 2.2 ] 
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[ Plot 2.3 ] 



32 



[ Plot 2.4 ] 




3. APPROACH FOR PATH AND FUNCTION 
GENERATION 


In this chapter we shall discuss a least square technique for synthesizing 4R mechanisms 
for path generation and function generation ( or coordinated input-output ) problems. This 
approach is generally used for more than three design positions. 

3 . 1 . PATH GENERATION 

Here the coordinates of some points on the specified path are available. These points may 
or may not be coordinated with the input angles. If the coordinated input angles are not 
specified then one can assume these angles which are required in this method. Now we 
have to synthesize a 4R mechanism such that its coupler point approximately passes 
through the specified points. 

Figure 3.1 shows a 4R mechanism, O2ABCO4 with unknown design parameters 
Ri, R 2 , R3, R4, R5, R6, a, p , </> and iff out of which a and j5 are assumed initially. 
Suppose there are n points through which the coupler point C has to pass and the coupler 
point at the i lh position is specified by its polar coordinates ( n, 8 i ) with the coordinated 
input represented by 021 - 

Now, writing the loop-closure equations for the loop OO 2 AC 1 O, we get , 

Ri cos a +R 2 cos 0 2 i + R 3 Cos y/\ = n cos 5; [3.1.1 ] 

and 

Ri sin a + R 2 sin 0 2 i + R 3 -sin iff i = n sin 8 \ . [3.1.2] 


Similarly, for the loop OO4BC1O, we have 


R 4 cos /3 +R 5 C 0 S 04 i + Recos i = r; cos <5 S 


[3.1.3] 
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and 


R 4 sin P +R 5 sin 0 4i + R 6 sin (j)\ = q sin <5 


Eliminating y/ ; from Equations 3.1.1 and 3.1.2, we get 


( Ri cos a + R 2 cos 6 2 i - n cos 8 ; f + ( Ri sin a + R 2 sin 6 2 ; - n sin 8 i 
or 


r 2 + R 2 2 + Ri 2 + 2 Ri R 2 cos ( 6 2 i - a ) -2 Ri n cos (a - 8 i ) 

- 2 R2 r, cos ( 0 2 i - 8 i ) = R3 2 


or 


R, [ 2 n cos (a - ) ] + R 2 [ 2 q cos ( 0 2 i - <5i ) ] + ( R 3 2 - R 2 2 - 

— Y\~ + R[ R? [ 2 COS ( 0 2i - cc ) ]. 


Similarly, eliminating 0 4 i from Equations 3.1.3 and 3.1.4, we get 

R 4 [ 2 r; cos ( P - 5 i )] + R 6 [2r i cos( 5 i )] + (R 5 2 -R 4 2 - 

= q 2 + R 4 R6 [ 2 cos ( 0 i - P )]. 


Substituting 


D, = Ri, D 5 = R 4 , 

D 2 =: R2> ^ ^6 = R-6* 


D 3 = R 3 2 - R 2 2 - Ri 2 , D 7 = R 5 2 - R 4 2 - Re 2 

D 4 = Ri R 2 , and D g = R 4 R 6 

Equations 3. 1 .5 and 3.1.6 are linearised in the design parameters Dj’s ( i 


D| [ 2 r, cos (a - ) ] + D2 [ 2 r, cos ( 02 - 5 i ) ] + D3 

- Tj 2 - D 4 [ 2 COS ( 0 2i - ® )] = 0 


[3.1.4] 


) 2 = R 3 2 


. 2 ) 

[3.1.5] 


h 

[3.1.6] 


[3.1.7] 

1 to 8 ) as 

[3.1.8] 


and 
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D 5 [ 2 q cos ( P - S i ) ] + D fi [ 2 q cos ( (pi - §i ) ] + D 7 

-rr-D 8 [2 cos ( (pi -/?)] = 0. [3.1.9] 

The design parameters are nothing but functions of the link lengths. If the coupler point 
passes exactly through the prescribed points at the desired input angles then Equations 
3. 1 .8 and 3. 1 .9 are perfectly satisfied. Here, only q, 5; and 0 2i (i = 1 ton) are available, 
but while solving Equations 3.1.8 and 3.1.9, it requires some assumptions 

( say a , ft and (p i ). At this stage, it must be appreciated that Equations 3.1.8 and 3.1.9 
are not satisfied exactly. Now let en and e 2 ; represent the error in the loop-closure 
equations 3.1.8 and 3.1.9, respectively, (at the i th design position ). Thus, 


eu = Di [ 2 q cos (a - ) ] + D 2 [ 2 q cos ( 0 2 i - <5, ) ] + D 3 

-rf-D 4 [2cos( 0 2i - a )] [3.1.10] 


and 


e 2 i = Ds [ 2 q cos (/? - §i ) ] +D 6 [ 2 q cos ( - <5* ) ] + D? 

-q 2 -D 8 [2cos( (pi- ft )]. [3.1.11] 


It is obvious from Equation 3.1.7 that 

D 4 = DiD 2 and D 8 = D 5 D6- 

Thus D 4 and D« make the Equations 3.1.10 and 3.1.11 non-linear. Let us first solve 
Equation 3.1.10 as explained below. 

Let 

D 4 = D,D 2 =X [3-1.12] 


and Dp = l p + X m p , for p — 1 , 2, 3. 


[3.1.13] 
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From Equations 3.1.12 and 3.1.13 

A =(\i+ X mi ) ( lo + A m 2 ). [3.1.14 ] 

Therefore, from equation 3.1.10 

on — ^ii,i -t" A eii, m , [3.1.15] 

where 

en,i = h [ 2 q cos (a - <5j ) ] + 1 2 [ 2 r, cos ( 02 i - 5i ) ] + 13 - n 2 

and 

en.m = mi [ 2 r ; cos (a - 8\ ) ] + m 2 [ 2 q cos ( 0 2i - Si ) ] + m 3 - 2 cos ( 6% - & )- 

In the ideal case, en = 0. To achieve this we would like to get l’s and m’s such that 
eiu = 0 and en, m = 0. 

Similarly, for solving Equation 3.1.11, 
let 

D 8 = D 5 D 6 = 77 
and 

D p = Up + 77 v p , for p = 5, 6, 7. 

From Equations 3.1.16 and 3.1.17, we get 

77 = ( ui + 77 vi ) (u 2 + 77 v 2 ). 

Therefore, from equation 3.1.11 

e 2 j = C 2 i,u + A e 2 j,v , 

j u = u, [2n cos (/3 - <5i)] + u 2 [2ri cos ( (pi- <5 i ) ] + u 3 - r; 


[3.1.16] 

[3.1.17] 

[3.1.18] 

[3.1.19] 


where 
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and 

e 2i,v = V, [ 2 H COS ( P - 5i)] + V 2 [2riCOS( fc - ) ] + v 3 - 2 cos ( <j>i-p ). 

Again in the ideal case e 2 j = 0. To achieve this we would like to get u’s and v’s such that 
eoj.u = 0 and e 2 i iV = 0. 

Using the least square technique, now we shall minimize the overall error ( e ), 
which is expressed as 

n 

e = X wh ( e H ,i 2 + e li>m 2 + e 2i , u 2 + e 2i , v 2 ), [3.1.20] 

(=1 


where the error functions are weighted with wtj for the i th position. 

It helps the designer to minimize the error more at the positions where he / she needs to 
do so. In other words, after some number of iterations, the designer can obtain a better 
solution with the help of these weighting functions wtj. In the first iteration, wt ; are taken 
as unity for all i’s and then it is checked at which design positions errors are more than 
acceptable. Thereafter, the wt;’s are modified accordingly. After the first iteration, we 
have taken wtj equal to the structural error at the i lh position. By iteration with these 
weighting functions the user can minimize the structural error at some critical points. 

Now, for the minimum value of e, we set the first partial derivatives to zero. 
Thus, we get 

Pj n pn 

— [ X wtj eiij 2 ] = — [ X wt;e lu 2 ] = 0, 
dh i= 1 d h i = 1 

’ [3.1.21a]' 

ph o d n i 

[X wtieu, in 2 ] '= - — [X Wti e H ,m~] = 0, 

dm 2 i=[ om 3 /=1 


P) n 

[ X wtj eii, 2 ] 


dl 


i = 1 


,m \ i = l 


[ X Wtj Cli.in - ] — 


[3.1.21b] 
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d n d n d n 

— - [ Z wtj e 2i ,u 2 ] = T~ [ Z wtj e 2i , u 2 ] = 3- [ Z wti e 2i , u 2 ] = 0 

du i i = 1 du 2 i= i du 3 i=1 


[3.1.21c] 


d n d n i d n 

— - [ Z wti e 2iiV 2 ] = 3- [ Z wti e 2i>v 2 ] = 3“ [ Z wti e 2i>v 2 ] = 0. 

dv\ /=1 dv 2 i=[ dv 3 i=1 


'3 i = i 


[ 3.1.21d ] 


Now, onwards the summation indices i and n are omitted. 

Equation 3. 1 .2 1 a gives 

1 1 Z [ wtj { 2riCos (a- 5i) } 2 ] 

+ 1 2 Z [ wti { 2 cos ( 0 2 j - <5i ) } { 2 n cos (a - <5 j ) } ] 
+ I 3 Z [ wti { 2 n cos {& - <5j ) } ] 


Z [ wti r 2 { 2 n cos (a - 8 \ ) } ], 


I, Z [ wtj { 2rjCos (a - 8\) } { 2r;cos( 9 2 \- S \) } ] 

+ I 2 Z [ wti { 2 r; cos ( 9 2 i - ) } 2 ] 

+ I 3 Z [ wt; { 2 Ij cos ( 02 i - <5i ) } 1 

= Z [ wt; r 2 { 2 r; cos ( 9 2 i - 8 ; ) } ] 


[3.1.22a] 


[3.1.22b] 


1 1 Z [ wtj { 2 r; cos ( a - 8 i ) } ] 

+ 1 2 Z [ wt; { 2 n cos ( 9 2 i - S j ) } ] 
+ I 3 Z [ wt i] 

= Z [ wt i r i 2 ]• 


[3.1.22c] 
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Thus, these three equations can be solved by using Cramer's rule to determine unknowns 
1 1 , lo and I 3 . 

In a similar manner, Equation 3. 1 .2 1 b gives 

mil [ wtj { 2 r, cos (a- <5; ) } 2 ] 

+ m 2 X [wti { 2 q cos ( 0 2i - Si) } { 2 r;cos(a- Si) } ] 

+ m 3 X [ wtj { 2 n cos (a - Si ) } ] 

= X [ wtj { 2 cos ( 0 2 i - a ) } { 2 q cos (a- ) } ], 

[3.1.23a] 

mi X [ wti { 2 r, cos (a- S t ) } { 2 r, cos ( 6 2 \ - S { ) } ] 

+ m 2 X [wti { 2 r iC os( 0 2 i - <5i) ) 2 ] 

+ m 3 x [ wti] { 2 rj cos ( 0 2 i - s i ) } 

= X [wtj [2 cos ( 0 2 i - & ) } { 2 Ti COS ( 02 i--Si) } ] 

[3.1.23b] 

and 

mi X [ wq { 2 q cos (a- 5; ) } ] 

+ m 2 X [ wti { 2 q cos ( 0 2 i - S i ) } ] 

+ m 3 X [ wh] 

= X [ wtj { 2 cos ( 02 i - & ) } ]• [3.1.23c] 

Again these three equations can be solved by using Cramer s rule to determine irq, m 2 
and m 3 . Once 1 1 , I 2 , I 3 and mi, m 2 , m 3 are known, we can determine X ( from Equation 
3.1.14 ), and then D t , D 2 and D 3 ( from Equation 3.1.13 ) and hence R|, R 2 and R 3 ( from 
Equation 3.1.7). 

Now we can determine y/ i’s, i=l to n, by using Equation 3.1.1 or 3.1.2. 

Assuming a suitable value for (j) 1 , we can determine (f>i ( for i = 2 to n ) by using 


the following relations. 
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0 i = 01 + (0i -!//•;), (for i = 2 to n ) [ 3.1.24] 

Thus, we have two values of A ( say A i and A 2 ) and correspondingly two values 
of Rt, R2 and R3 ( say Ri[l], R 2 [l], R3[l] and Ri[2], R 2 [2], R3[2] ) and hence two values 
of 0i’s ( for i = 2 to n ). 

In the same manner, Equations 3.1.21c and 3. 1.2 Id give 


u, I [ wti { 2 q cos (/3 - 8 \ ) } 2 ] 

+ u 2 X [ wq { 2 q cos ( 0 ; - 8 i ) } { 2 q cos ( /3 - <5 i ) } ] 

+ u 3 X [ wtj { 2 q cos ( ft - 8 i ) } ] 

= X [ wtj q 2 { 2 q cos (/3 - 5i ) } ], 

[3.1.25a] 

ui X [ wti { 2 q cos (J3 - 5i ) } { 2q cos ( 0 i - $i ) } ] 

+ u 2 X [ wtj { 2 q cos ( 0 •, - 8 j ) } 2 ] 

+ u 3 X [wti { 2qcos( 0i- 5i) } ] 

= X [ wt; q 2 { 2 q cos ( <p i - 8i ) } ], 

[3.1.25b] 


Ui X [ wt; { 2qcos(j8 - 5j) } ] 

+ u 2 X [ wtj { 2 q cos ( 0i - 8i ) }] 
+ u 3 X [ wt; ] 

= X [wt; q 2 ], 


[3.1.25c] 
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vi X [wtj { 2 n cos (/3 - Si) } 2 ] 

+ v 2 X [wtj { 2r;cos( (pi - S t ) } { 2rjCos(/f- <50 } ] 

+ v 3 X [ wtj { 2 n COS ( P - S i ) } ] 

= X [ wtj { 2ijCos( 0i- p ) } { 2r,cos(jB-Si) } ], 

[3.1.26a] 


and 


V, X [ wt| { 2 n cos C/3 - > } { 2 n cos 

+ v 2 X [wti { 2ncos( 0i- Si) } 2 ] 

+ v 3 X [ wtj { 2 n cos ( - Si ) } ] 

= X [ wtj { 2 r; cos ( 0; - p ) } { 2 rj cos ( 0; - Si ) } ] 

[3.1.26b] 


v, X [wt; { 2 rj cos (/3- Sj) } ] 

+ v 2 X [ wtj { 2 n cos ( (j> i - S i ) } ] 

+ v 3 X [ Wtj ] 

= X [ wti { 2 n cos ( 0 i - >5 ) } ]. 


[3.1.26c] 


Again, from Equations 3.1.25a, 3.1.25b and 3.1.25c we determine ui, u 2 and u 3 . 
Similarly Vi, v 2 and v 3 are obtained from Equations 3.1.26a, 3.1.26b and 3.1.26c. 
Therefore we get r\ ( from Equation 3.1.18 ), and then D4, D5 and D6 ( from Equation 
3.1. 17 ) and hence R4, R5 and R 6 ( from Equation 3.1.7 ). There are two values of <p\S ( 
for i = 2 to n ) and each value of 0 ;’s ( for i = 2 to n ) gives two values 77 [ say the first 
value of 0j’s ( for i = 2 to n ) gives 77 n, 77 12 and the second value of 0i’s ( for i = 2 to n 
) gives 77 21 , 77 22 ]. Hence there are four values 77 and therefore, we get four values of R4, 
R5 and Rf„ which are denoted as 
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R4[1][1],R4[1][2],R4[2][1],R 4 [2][2], 
R5[l][l], Rs[1][2], R S [ 2][1 ], R s [2][2] 

and 

R6[1][1],R 6 [1][2],R 6 [2][1],R 6 [2][2]. 
Hence, we have four solutions as 


[i] R 1 [l],R 2 [l],R 3 [l],R 4 [l][l],R 5 [l][l] ) R 6 [l][l] a, p, <f> u 

[ii] R l [l],R 2 [l] 1 R 3 [l],R4[l][2],R 5 [l][2],R 6 [l][2] a, p, <p { , 

[iii] Ri[2], R 2 [2], R 3 [2], R 4 [2][l], R 5 [2][l], R 6 [2][l] d, p, 0, 

[iv] R[[2], R 2 [2], R 3 [2], R 4 [2][2], R s [2][2], R 6 [2][2] dj , 0,. 


Now, the parameters ( R|, R 2 , R 3 , R 4 , R 3 , R a , P and 0i ) can be converted 
into the link lengths ( L|, L 2 , L 3 and L 4 ) and the polar coordinates ( R, y ) of the coupler 
point, such that 

Li = 0 2 0 4 = [ ( R 4 cos P - R| cos a ) 2 + ( R 4 sin P - Ri sin a ) 2 ] 1/2 , 

L 2 = o 2 a = R 2 , 

L 3 = AB = [ R 6 2 + R 3 2 - 2 R fi R 3 cos ( 0 , - \jf , ) ] 1/2 , 


l 4 = o 4 b= R 5, 


[3.1.27] 


R = AC = R 3 , 

and 7 = Z CAB = sin"‘[ R6 sin ( $ \ - 17 1 ) / L 3 ]. 

These new parameters ( L,, L 2 , L 3 , L 4 , R and 7 ) are used in the animation of a 4R 
mechanism for path generation. 
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Hence, we can find L,[k], L 2 [k], L 3 [k], L 4 [k], R[k] andy [k] for each k ,h solution. 
It is not guaranteed that all the above solutions give acceptable results. It mainly depends 
upon the assumed values of a , /? and tp , needed to solve the equations. By changing 
these values one can arrive at better solutions. 

Flow Chart - Chart 3.1 shows the steps used in the programming. Here are some 
mathematical bridges ( Quadl, Sqrtl, Quad2, Sqrt2, see in Chart 3.1 ), without getting 
these positive one can not move further. These expressions are controlled by changing the 

assumed values of &, p and <p \. Some bridges are being controlled automatically and 

others need manual controlling. Regarding the selection of a and ft , one should choose 
the parameters such that most of the specified points lie near the perpendicular to the line 
O2O4 and to get a crank-rocker or a double-rocker mechanism, the specified points must 
lie on one side of the line O2O4. 

The present software can also handle some constraints like 

[a] the ratio ( L max / Lmin ), where L max and L min are the maximum and minimum values of 
link lengths of the mechanism, 

[b] the designed mechanism must be a crank-rocker and 

[c] the structural error at a specified point must be zero i.e, the coupler point must pass 
through a specific point on the path. 

The ratio L max / L min is controlled by j8 whereas a crank-rocker can be ensured by 
(j) |. Here, to minimize the structural error, the least square error function is weighted by 
the distance between the specified and generated points at each prescribed 9 2 . Obviously 
some trials are necessary to arrive at an acceptable solution. Here some steps are given, 
following which one can easily run the program quickly to get a good result. 


[ i ] Give n, (Xj, y-, ) or ( n, S ■, ) and 6 r,. 
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[ ii ] Select any one of the solutions ( 1, 2, 3 and 4 ) and assume {a, ft and Q x •) as 
default values. Initially the constraints are left out. 

[ iii ] Quad 1 is automatically controlled, but controlling of Sqrtl may require a manual 
change of a. Similarly, Quadl is controlled automatically but Sqrt2 is controlled by 

changing p manually. The program itself takes care of the constraints on the ratio ( / 
L in i n ) and the crank-rocker compulsion. 

[ iv ] Check the other three solutions and choose the best one. Now, if the structural errors 
: are large, then go to step [ii] and change the assumptions keeping in mind the previous 
results and the assumptions made. After two to three trials, one can guess the suitable 

values of a , p and <p The variation of 0 ( can help to reach a good solution. 

[ v ] For better results, choose an option of weighting the error function. Weighting on the 
error functions are possible only in the case of crank-rocker and double-rocker 
mechanisms. One should check the other three solutions if one solution exists. 

If instead of prescribed points, a path is given, then start solving the problem by 
considering 6 points on the path at equal steps of 6o- Then increase the number of points 
to 8 , 10 and 12. The present software considers a maximum of twelve points but one can 
replace it by any number. However, care should be taken regarding the memory and the 
speed of the computing system. 

An Illustration - Let us consider a path generated by a crank-rocker mechanism ( shown 
in Plot 3.1 ) with 

Li = 2 unit, 

L 2 =1 unit, 

L 2 = 2 unit, 

L 4 = 2 unit, 

R = 3 unit 
7 = 120°. 


and 



Number of points selected = 

: 12, which are 



Point 

X 

coordinate 

y 

coordinate 


.8906 


-2.0701 


-2.2042 


-2.5835 


-3.0845 


-3.5495 


-3.8435 


-3.8451 


-3.4730 


-2.7701 


-1.9771 


-1.5689 


.8028 


1.1157 


2.1651 


2.525 


2.3894 


1.8414 


0.9564 


-0.1465 


-1.2682 


-2.1518 


-2.5584 


-2.2525 


0.0 


30.0 


60.0 


90.0 


120.0 


150.0 


180.0 


210.0 


240.0 


270.0 


300.0 


330.0 


For this problem, assuming a , j3 , (j) \ and L max / L m i n as default values (a = 56°, ft 
354°, (j) i = l 66° and L max / L min =20 ) we get the following results. 


Solution 

no. 


a 

A 

P 

0i 

Linax 

! L m in 

Type of 

mechanism 

( NP )Number of 

specified points through 

which the coupler point 

approximately passes. 

56 

354 

166 

19.62 

Crank - rocker 

2 

56 

354 

166 

1.497 

Non Grashofian 

6 

56 

354 

166 

10.13 

Non Grashofian 

2 

56 

354 

166 

2.289 

Non Grashofian 

8 
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In the next iteiation, an attempt is made to get a crank-rocker mechanism for 
solution no. 4. The final results turn out as L IIiax / L min = 1.968 and NP = 4 for a crank- 

rocker mechanism. Then for a new set of assumptions (a = 50°, /? = 50°, <p , = 90° and 
L ma x / Lmin =20 ), the following results are obtained. 


Solution 

no. 

a 

P 

0i 

max 

/ L min 

Type of 

mechanism 

(NP)Number . of 

specified points through 

which the coupler point 

approximately passes. 

i 

50 

50 

90 

13.56 

Crank - rocker 

5 

2 

50 

50 

90 

1.852 

Non Grashofian 

3 

3 

50 

50 

90 

13.07 

Crank - rocker 

4 

4 

50 

50 

90 

5.366 

Crank - rocker 

9 


Out of all the four solutions, the fourth one can be selected ( see Plot 3.2 ). To 
observe the accuracy of the results, X_error ( % ) and Y_error ( % ) are introduced. These 
errors, at an i lh position ( i.e. corresponding to 6o\ ), are defined as 

( % ) X_error [i] = [ { ( Xj ) g - ( Xj ) s } 100 / { ( x max )s ■ ( x m m )s } ] 
and 

( % ) Y_error [i] = [ { ( yi ) g - ( Yi )s } 100 / { ( y m ax )s - ( Yinm ) s } ]» 


where i = 1 to n ( no of specified points ), 

suffixes ‘g’ and ‘s’ denote the generated and specified points, 

Xmax = the maximum value of Xi ( x _ coordinate of the point at 

the i th position ), 

x m in = the minimum value of Xj, 

ymax = the maximum value of y, ( y _ coordinate of the point at the i 
position ) 
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and y,nin = the minimum value of yj. 

To reduce the errors near the points 3 and 4 ( as expected ) weighting of the error 
function has been considered and after the first iteration a better solution (regarding those 
two points ) is obtained ( see Plot 3.3 ). The errors at all the points corresponding to 0 2i ( 
i = 1 to n ), before and after weighting the error functions, are given below. 



Before weighting the error 

function 

After weighting the error 

function 

Point 

e 2 

( % ) X_error 

( % ) Y_error 

( % ) X.error 

( % ) Y_error 

1 


14.82 

-15.19 

6.22 

-7.83 

2 

30 

-2.65 

-14.99 

-2.50 

-8.11 

*3 

60 

21.41 

14.21 

14.11 

10.30 

*4 

90 

13.93 

8.67 

7.76 | 

5.23 

5 

120 

3.029 

-2.19 

-1.24 

-0.96 

6 

150 

-2.70 

-2.58 

-5.27 

-5.38 

7 

180 

-3.07 

-4.73 

-4.20 

-6.99 

8 

210 

-0.71 

-3.62 

-0.86 

-5.10 

9 

240 

-0.0002 

0.0065 

0.14 

0.43 

10 

270 

-3.28 

3.925 

-3.92 

4.79 

11 

300 

1 

-7.34 

5.414 

-10.1 

7.85 

12 

300 

-1.63 

1.039 

7.83 

5.485. 


The errors at the other points ( e.g. 1 1 and 12 ) are now more. So applying weights 
is more useful when the structural error at some points are critical for a particular 


application. 
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3.2 . FUNCTION GENERATION 

If a 4R mechanism has to generate a particular function (which is given ), with more than 
three precision positions, then the optimization method is more useful regarding the 
overall minimization of the error. 

Let 9 2-oi and 9 4 -oi be the input and output angles, respectively at the i th position as 
shown in Fig 3.2. Now let us assume the initial positions of the input and output lin- 
ks, which are at angles 9 2-0 and 9 4-0, respectively ( Fig 3.2 ). Therefore, from Fig 3.2, for 
n design positions 

9oj= 02-0+ 9 2 - 0 i ( i = 1 to n ) [3.2.1] 

and 9 4 j = $ 4 -o + 0 4 -oi ( i = 1 to n ). [ 3.2.2 ] 

Let 0 4 i be the actual angular displacement of the output link of the synthesized mecha- 
nism corresponding to the input angle 9 2 i- Writing the well known Freudenstein’s 
equation m for the i th position, we get 



D, 

COS 6 2\ ~ D2 COS 0 4i 4- D3 — COS ( 9 2i 

" 04i )j 

[ 3 . 2,3 ] 

where 






D, 

= Li / L 4 , D 2 = L,/L 2 



and 




[3.2.4] 


D, 

= ( L|“ + L 2 2 - L3 2 + L4 2 ) / ( 2 L 2 L4 ). 




If the angle 0 4i is replaced by 9 A \ ( which is available at 9u ) then Equation 3.2.3 is not 
satisfied in general and the difference between the L.H.S. and the R.H.S. of Equation 
3.2.3 gives a small error. While using the least square technique, let e be the error 

defined as 

/ 

e = X [ Di cos 0 21 - D 2 cos 9 a\ + D 3 - cos ( 9 2 \ - 04 i ) f [ 3 - 2 - 5 3 

i = l 
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Now, for the error e to be minimum, we take the partial derivatives of e with 
respect to D|, D 2 and D 3 separately, and set each of the derivatives to zero. Thus 


de 

<9Z), 


n 

2 X [ Dj cos 0 2i - Do cos 9 4 j + D 3 - cos ( 9 2l - 9 4 ; ) ] cos 9 2 i 
i=\ 


= 0, [ 3.2:6 ] 


de 

Jd~ 2 


n 

-2 X [ Di cos 9 2i - D2 cos 0 41 + D3 - cos ( 6 n - 0 4i ) ] cos 0 4 i 
i = 1 


=0 


[ 3.2.7 ] 


and 


d e 


n 

2 I [ Di cos 0 2 i - D 2 cos 9 41 + D 3 - cos ( 0 2 i - 9 4 i ) ] 

i'=l 


= 0 . 


[3.2.8] 


The above equations are rearranged to obtain 


n n n 

D| S COS 2 9 2i - D 2 X COS 0 4i COS 02i + ©3 Z COS 02i 
i — 1 i= 1 1=1 

= X cos ( 02i - 04 i ) COS 02i, [3.2.9] 

i= 1 

n n n 

Di X COS 0 2i COS 0 4 i - D 2 X COS" 0 4 i + D 3 X COS 0 4 i 
I=1 i=l 1=1 

= X COS ( 02i - 04i ) COS 04i [3.2.10] 

j=l 


n n ^ 

and Di X cos 0 2 i - D 2 X cos 0 4 i + D 3 X 1 
i=l i= 1 i=1 

= X COS ( 9 2\ ~ 0 4i ) [3.2.11] 

f=l 

The above three equations are linear, simultaneous and non homogeneous in three 
unknowns, Di, D 2 and D 3 . Hence, these can be solved by using the Cramer s rule. 
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Finally from Equation 3.2.4 we get the link lengths. By changing the assumed 
values of 02.o and ««, we can reduce the absolute value of ( f>„ - 0 4i ) a t each design 
position and thus get a better solution. 

Flow chart with an illustration - The programming for the optimized function 
generation problem can be easily understood from Chart 3.2. There are mainly two 
variables 

( 0 2 -o and 0 4 -o ) to be assumed and one more parameter Li is used to define a unique 
mechanism. An induced condition ( Sqrt3 > 0 ), mentioned in the Chart 3.2, can easily be 
satisfied. For obtaining a better solution, one should minimize I 6 4 . 0 - 0 4i I by changing 
the initially assumed values. After computing the link lengths, DRAW_4R( ) is executed 
to animate the mechanism. For verification of the result we print I 0 4i - 0 4i I at each 0 2 > 
and one can also display the result graphically. 

We have considered a problem for which a set of 8 coordinated input-output 
angles are given. These are given below. The angles are expressed in degrees. 


i 

!■ 

2 

3 

4 

5 

6 

7 

8 

0 2-0i 

15 

60 

105 

150 

195 

240 

285 

330 

0 4-0i 

30 

50 

70 

90 

110 

100 

85 

75 


Assuming 0 2 -o = 0.0, 0 4 -o = 0.0 and Li = 5, the dimensions of the synthesized 
mechanism along with the plot of On vs. (j >4 are shown in Plot 3.4. 

The errors ( i.e., I 0 4 j - 0 4 i I ) at different positions are given mentioned below. 


i 

n 

2 

3 

4 

5 

6 

7 

8 

0 2-0i 

15 

60 



195 

240 

285 

330 

1 0 4i - 0 4i 1 

■ 

2 

0 

0 

6 

6 

10 

12 
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[ OPTIMIZATION APPROACH FOR PATH GENERATION ] 
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Change f) 


'Quad2*>0 


•compute 77 [1], 77 [2] from Equation 3.1.18 
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[Chart 3.1 ] 











Compute D| , D 2 & D 3 from 
Eqns 3.2.9, 3.2.10 & 3.2.11 



[ Chart 3.2 ] 
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[Plot 3.1 ] 
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[Plot 3.2] 
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[Plot 3.3 ] 
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[ Plot 3.4 ] 




4. COUPLER CURVES 


In this chapter, two important aspects of 4R coupler curves are presented. The first one is 
related to the generation of two cusps by a coupler point of a4R mechanism. The second 
one is concerned with the coupler curve cognate linkages. 

4.1 TWO CUSPS PROBLEM 

The objective is to design a crank-rocker mechanism, in which a coupler poirit has two 
cusps at prescribed locations and the amount of crank rotation as the coupler point moves 
from one cusp to the another is also prescribed. 

Before solving this problem, we should note the following result [1] . 

For every crank-rocker mechanism, there exists a coupler point which has two 
cusps. This coupler point can be located at the instantaneous center corresponding to the 
configurations when the crank-pin A lies on the Thales circle CR T ( a circle with O 2 O 4 as 
a diameter) (Fig 4.1 ). 

It can be easily proved that the sum of the crank rotation (6 2c ) and the rocker 
rotation ( 6 4c ) as the coupler point goes from one cusp to the another is % . This can be 
proved from Fig 4.1 as explained below. 

Since A 0 4 B 2 A 2 is obtained by rotating A 0 4 BiAi about 0 4 through an angle 
0 4c and Z B 2 0 4 Bi = Z C 2 0 4 Ci = 6 4c , we get Z A 2 0 4 Ai = 0 4c . 

But A|, A 2 , 0 2 and 0 4 are on the same circle CR T , therefore 

02c+04c=7r or 64 C =n-9 2c . [4.1.1] 

Let Ci = (xi, yi ) and C 2 = (x 2 , y 2 ) are the prescribed locations of the cusps 
whereas 0 2c is the given crank rotation as the coupler point goes from one cusp to the 
another. 

It is also known, that the double point ( if it exists on the coupler curve ) lies on a circle 
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known as circle of singular foci 1 ' 1 . But, the circle of singular foci passes through the 
fixed hinges ( i.e. 0 2 and 0 4 ). Hence we can say that Q, C 2 , 0 2 and 0 4 will lie on the 
same circle. 

From Fig 4. 1 , concentrating on A ’s 0 4 AiCi and 0 4 A 2 C 2 , we observe that 


Z Ai0 4 C| — ZA|0 4 B 2 + 6 4 C = Z A 2 0 4 C 2 , 


Z CiAi0 4 = Z C 2 A 2 0 4 = K / 2 
and 

0 4 Ai = 0 4 A 2 . 

Therefore, we can say that A 0 4 AiCi = A 0 4 A 2 C 2 and 0 4 Ci = 0 4 C 2 . 

Thus, for the synthesis of the desired mechanism, we draw a circle ( CRi ) with 
the following characteristics : 

[i] It passes through Ci and C 2 . 

[ii] The angle subtended by C)C 2 at 0 4 = ( K - 0 2c ). 

[iii] 0 4 C, = 0 4 C 2 . 

In other words, we can draw a circle CRi passing through Q, C 2 and 0 4 such that 0 4 is 
on the mid-normal of CiC 2 and the angle subtended by C]C 2 at 0 4 is ( 7t - 62 c ) ( Fig 4.2 

)• 

Suppose 

Co = ( xco, yco ) is the center of CR t ( of radius ri ), 

0 4 = ( X04, y04 )> 
p = the distance of Co from CiC 2 
and Ci 2 = distance between Ci and C 2 . 

Therefore, from Fig 4.2 

C\ 2 = sT [ ( xi - x 2 ) 2 + ( yi - Y 2 ) 2 1, [4-1.2] 


8 ( slope of CiC 2 ) = tan' 1 [ ( yi - y2 ) / ( Xi - x 2 ) ], 


[4.1.3] 
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i' i = C 12 / ( 2 sin 0 2 c X 

[4.1.4] 

P=V~ C r i 2 " ( C| 2 2 / 4 ) ], 

[4.1.5] 

Xco = X 2 + (Ci 2 / 2 ) cos S + p sin 8 , 

[4.1.6a] 

yco = y 2 + ( C 12 / 2 ) sin ^ - p cos 8, 

[4.1.6b] ■ 

xo4 = X 2 + ( C [2 / 2 ) cos 8 + ( p + ri ) sin <5 

[4.1.7a] 

and 


yo4 = y 2 + ( C 12 / 2 ) sin 5 - ( p + ri ) cos 8 . 

[4.1.7b] 

Writing the equation of the circle CR|, we have 


( x- xco f + ( y- yco ) 2 = r, 2 . 

[4.1.8] 

Since 0 2 [ = ( xq 2 , yo 2 ) ] also lies on the same circle CRi, 

X 02 can be chosen 

anywhere on the circle between ( xco - H ) and ( xco + ri ). Hence, from Equation 4.1.8 

yo 2 is known. Thus O2 and O4 are known, and we can draw the Thales circle CR T . 

Let Ct [ = ( Xct, ycr ) ] and r t be the center and radius of the circle CRt, 

respectively. 

Therefore, xct = ( X 02 + X 04 ) / 2j 

[4.1.9a] 

ycT = ( yo 2 + yo4 ) / 2 

[4.1.9b] 

and 


r t = [ ■yj~ { ( X 04 - X 02 )” + ( yo 4 - yo 2 ) 2 } ] / 2. 

[4.1.10] 


Equation of the circle CRt is given by 
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( x- x C T f + ( y- y CT ) 2 = r, 2 . 


The equation of the line O2C1 is 
y = m x + c, 

where 

m = (yi -y 02 )/(xi -x 02 ) 
and 

c = y 0 2 - m X02 


[4.1.11] 


[4.1.12] 


[4.1.13a] 


[4.1.13b] 


We know that the crank-pin lies on the point of intersection of the crank circle and 
the Thales circle when the coupler point lies at either of the cusps. A cusp is an 
instantaneous center of velocity of the coupler. Also for a 4R mechanism, the 
instantaneous center of the coupler relative to the fixed link O2O4 lies at the intersection 
of lines drawn through the input and the output links ( using the well-known theorem of 
three centers ). In other words, at the configuration when the coupler point lies at the 
instantaneous center ( or cusp on the coupler curve ), the crank-pin lies on the line joining 
the fixed end of the crank and the instantaneous center ( or cusp on the coupler curve). On 
this basis, we can say that the crank-pins Ai [ say (x A [, yAi ) ] and A 2 [ say (x A2 , ytci ) ] ( 
lying on the Thales circle CRt ) lie on the lines 0 2 Ci and 0 2 C 2 respectively ( Fig 4.2 ). 

Now, from Equations 4.1.11 and 4. 1 . 1 2 


( x- Xct ) 2 + (m x + c - ycr ) 2 = r t 2 
or 

( 1 + m 2 ) x 2 + 2 ( ti + m t 2 ) x + ( ti 2 + 1 2 " -r t 2 ) = 0, [ 4.1.14 ] 


where 


t, = - xcr and t 2 = c - ycr- 


The quadratic Equation 4.1.14 gives two values of x, as the line O2C1 cuts the circle CRt 
at two points. Solving this equation, we get the x coordinate of Ai i.e. x A i- Then, from 
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Equations 4. 1 . 12, 4. 1 .13a and 4.1.1 3b we can find the y _ coordinate of Ai i.e. y^i. Hence 
A| is known. 

To find all the link lengths, we have to locate the point Bi [ say ( x B i, y B i ) ]- 
Since the point Bi (the free end of the follower when the coupler point lies at the cusp Ci 
) lies anywhere on the line O4C1, Therefore, we can chose x B i between x ( ( the x _ 
coordinate of Q ) and X04 (the x _ coordinate of O4 ). 

Now the equation of line O4C1 is written as 

y = m ( x + ci, [4.1.15] 

where 

mi =(yi -yo 4 )/(xj -x 04 ) [4.1.16a] 

and 

ci =y 04 - mi x 04 . [4.1.16b] 

Therefore, we can find y B i corresponding to x B i on the line O4C1 ( from Equation 
4. 1 . 15 ). Depending upon the selection of x 0 2 and x Bi we can find many solutions but we 
control these selections to get a crank-rocker mechanism with 0 2 Ai as the shortest link. 
Thus a crank-rocker mechanism O2A1B1O4 is designed. 

Flow chart with an illustration - Chart 4.1 clearly explains the steps used to design a 
crank-rocker mechanism for the present problem. The assumptions made are 

X02 = xco - n /2 
and x B i = (xi + Xo 4 )/ 2 , 

which can be changed to get a required mechanism. As an illustration, we have selected 
two cusps at 

Ci = ( 10, 8 ), 

C 2 = (1,1) 
and 0 2c = 1 20°. 

The designed mechanism with coupler curve and link lengths are shown in Plot 4.1. 
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4.2 COGNATE LINKAGES 

The term cognate ( coined by Hartenberg and Denavit ) of a mechanism means a different 
mechanism which generates the same coupler curve. The most important theorem ni , 
developed by Samuel Roberts ( 1875 ) and Chebyschev ( 1878 ) indepen-dently, states 
that three different planar 4R linkages trace identical coupler curves. 

Let O 2 ABO 4 is a given crank-rocker mechanism and C is its coupler point ( Fig 
4.3). There exist two 4R cognate linkages of this mechanism. To get the cognate linkages, 
first draw parallelograms 0 2 ACD and O4BCG. On sides DC and CG, make triangles 
EDC and FCG respectively similar to the triangle CAB. Again draw a parallelogram 
OECF and make the hinge O fixed. This ten-link mechanism is found 1 - 11 constrained ( 
with degree of freedom = 1 ). Thus, the common coupler point C of mechanisms 
O 2 ABO 4 , O 2 DEO and OFGO 4 generates the same coupler curve. The mechanism 0 2 DE0 
is called the left-hand cognate and OFGO4 the right-hand cognate. 

It should be noted [11 that, A OO 2 O 4 is also similar to A CAB. 

Now, suppose 

O 2 O 4 = Li, 0?A = L 2 , AB = L 3 , 

0 4 B = L 4 , AC = R and BC = S. 

The unknown design parameters of both cognates can be found analytically in the 
following manner. 

[i] For OFCGO4 ( Right-hand cognate ) - The design parameters RG_Li, RGJU, 
RG_L 3 , RG_L 4 , RG_R and RG_S are given as below. 

From Fig 4.3, 

RGJLi = 00 4 = ( 0 2 0 4 ) ( BC/ AB ) 

( from similar A ’s OO2O4 and CAB ), 

or RGJLi = Li ( S / L 3 ), [4.2.1] 


RG_L 2 = OF = EC = LG_S = ( S / L3 ) L 2 , 


% 


[ 4.2.2 ] 
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RG JL 3 = FG = ( sin 7 1 / sin 73 ) GC 

( from A CDE ) 

or 

RG_L 3 = ( S / L3 ) L 4 

( from A ABC ), 



[4.2.3] 


rg_l 4 = 0 4 G = S, 

[4.2.4] 


RG_R = FC = ( sin 7 2 / sin 73 ) GC 

( from A CDE ), 

or 

RG_R = (R/L3)L 4 

( from A ABC ), 



[4.2.5] 

and 

RG_S = GC = L 4 . 

[4.2.6] 


[ii] For O2DCEO ( Left-hand cognate ) - The design parameters LGJU, LGJL 2 , LGJL 3 , 
LG_L 4 , LG_R and LG_S are given as below. 

From Fig 4.3, 

LGJU = 0 2 0 = (O 2 O 4 ) (AC / AB ) 

( from similar A ’s OO2O4 and CAB ), 


or 

LGJL, = Li ( R/L 3 ), 

[4.2.7] 


LG_L 2 = 0 2 D = AC = R, 

[4.2.8] 


LG_L 3 = DE = ( sin 72 / sin 7 3 ) DC 

( from A CDE ) 

or 

LG_L 3 = ( R / L3 ) L 2 

( from A ABC ), 



[4.2.9] 


LG_L 4 = OE = FC = RG_R = ( R / L3 ) L 4 , 

[4.2.10] 


LG_R = DC = OoA = L 2 

[4.2.11] 
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and 

LG_S = EC = ( sin 7 i / sin y 3 ) DC ( from A CDE ) 

°r LG_S = ( S / L3 ) L2 ( from A ABC ). 

[ 4.2.12 ] 

One can easily prove that the maximum-minimum link length ratio ( / L^,, ) 

for the original and its cognate linkages are same. 

We can analyze the angular relationships of the various links of different cognate 
linkages in the following manner. 

Observing the parallelograms OiACD, O4BCG and OECF we can say that 


Original Linkage 
Angular motion of link 2 = 
Angular motion of link 3 = 
Angular motion of link 4 = 


Left - Hand Cognate 
Angular motion of link 6 = 
Angular motion of link 5 = 
Angular motion of link 7 = 


Right - Hand Cognate 
Angular motion of link 8 
Angular motion of link 10 
Angular motion of link 9. 


Considering the above relationships we can conclude that, if the original linkage 
( O2ABO4 ) is a crank-rocker with say 2, as the crank ( when 3 and 4 must oscillate ), then 
the right-hand cognate is also a crank-rocker, with 8 as the crank ; however, the left-hand 
cognate then must be a double-rocker ( with the coupler 6 undergoing the full rotation ). 
Similarly both the cognates of a double-rocker linkage ( say 0 2 DCE 0 ) are crank-rockers 
( OFCGO4 and O2ACBO4 ). 


It should be noted that the transmission angle ( /x min ) for crank-rocker cognate 
linkages are same. 

From Equation 2.2. 1 , the transmission angle for the mechanism is 


/j, — cos' 1 [ { ( L3 2 + L4 2 ) - ( Li 2 + L 2 2 ) - 2 Li L 2 cos 6 2 } / ( 2 L3 L4 ) ]. 

[4.2.13] 
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The transmission angle for the right-hand cognate is given by ( replacing L’s by RG_L’s ) 

M right = cos' 1 [ { ( RG_L 3 2 + RG_L 4 2 ) - (RG_L, 2 + RG_L 2 2 ) 

- 2 RG_L, RG_L 2 cos 0 2r } / (2 RG_L 3 RG_L 4 ) ] 
where 6 2r represents the rotation of link 8. 

Substituting the values of RG_L,, RG_L 2 , RG_L 3 and RG_L 4 in the above equation 
( from Equations 4.2. 1 - 4.2.4 ), we get 

// right = cos' 1 [ [ { (L 4 S/L 3 ) 2 + S 2 } - { (L^S/LO 2 + (L 2 S/L 3 ) 2 } 

- 2 (L, S/L 3 ) (L 2 S/L 3 ) cos e Zl ] / { (L 4 S/L 3 ) S } ] 

or 

jj right = COS ' 1 [ { ( L 3 2 + L 4 2 ) - ( L, 2 + L 2 2 ) - 2 L, L 2 cos 0 2r } / ( 2 L 3 L 4 ) ]. 

[4.2.14] 

Suppose the original mechanism ( 0 2 AB0 4 ) is a crank-rocker. Since, link 8 is equiva- 
lent to link 2 regarding the angular motion, we can say that # 2r and 6 2 posses the same 
values for calculation of the minimum transmission angle. Hence, from Equations 4.2.13 
and 4.2.14, fright and fi are same. 


Illustration - For every crank-rocker and Grashofian double-rocker mechanisms, we can 
visualize the other two cognates using the present software. While solving a path genera- 
tion problem for three precision points, we get a double-rocker mechanism, but as per the 
requirement of a crank-rocker mechanism one can select one of its cognates. Plot 4.2 
shows the original mechanism ( double-rocker ) and its cognate linkages ( crank-rockers ) 
along with the coupler curve passing through the specified points in a path generation 
problem. 







o 



[Fig 4.3] 
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[ TWO CUSPS PROBLEMS ] 



[ Chart 4. 1 ] 
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[Plot 4.1 ] 
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5. KINEMATIC CONDITIONS 


In this chapter, we are going to present the synthesis of 4R mechanisms for prescribed 
kinematic parameters ( e.g. angular velocities and angular accelerations etc. ) of some of 
its links. Let 6 \, ( 0 { and a \ are the angular displacement, angular velocity and angular 
acceleration of the i th link as shown in Fig 5.1. Depending upon the values of i we can 
have two cases. 

Case 1. Tf 0 j, j and a \ for input ( i = 2 ) and output ( i = 4 ) links are given. 

Using Freudenstein ’s equation we have 
D| cos 6 2 + D 2 ( -cos 04 ) + D 3 = cos ( 02- 04 ), 

where 

Di = L| / L4, 

D 2 = L, / L 2 

and D 3 = ( L, 2 + L 2 2 - L 3 2 + L 4 2 ) / ( 2 L 2 L 4 ) 

with L as the length of the i th link, i = 1 to 4. 

Differentiating Equation 5. 1 with respect to time, we get 

Di ( (U 2 sin 0 2 ) + D 2 (- <S) 4 sin 0 4 ) = ( co 2 - CO 4 ) {sin ( 0 2 - 0 4 ) } 

[5.3] 

and differentiating Equation 5.3 once more with respect to time, we obtain 

Di ( co 2 cos 0 2 +a 2 sin 0 2 ) + D 2 (- tu 4 2 cos 0 4 - ausin 0 4 ) 

= ( co 2 - ® 4 ) 2 { cos( 0 2 - 0 4 ) } + ( a 2 - a 4 ) { sin ( 0 2 - 0 4 ) }• 

[5.4] 

Three design parameters, Di, D 2 and D 3 can be solved from Equations 5.1, 5.3, and 5.4 
using Cramer’s rule. Then using Equation 5.2, and assuming any one of the link-lengths ( 
i.e., the scale of the link diagram ) we can have the other three link-lengths. 


[5.1] 


[5.2] 
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Case 2. If co \ and a j for input ( i = 2 ), coupler ( i = 3 ) and output ( i = 4 ) links 
are given. 

Writing the horizontal and vertical components of the loop closure equation for 
the loop O 2 ABO 4 O 2 , we get 

L 2 cos 9 2 + L 3 cos 9 3 = Li + L 4 cos 9 4 [ 5.5 ] 

and 

L 2 sin 9 2 + L 3 sin 9 3 = L 4 sin 9 4. [ 5.6 ] 

Equations 5.5 and 5.6 can be combined into a complex exponential form as given below. 

L2 exp ( i 9 2 ) + L3 exp ( i 9 3 ) = Li + L4 exp ( i 9 a ) 
or 

( L 2 /Li ) exp (i 0 2 ) + ( L3 / Li ) exp (i 9 3 ) - ( L4 / Li) exp ( i 0 4 ) = 1 
or 

Z 2 + Z 3 - Z 4 = Z, [ 5-7 ] 

where Zn [ = (L„ /L, ) exp( i 9 n ), n = 1, 2, 3, 4 ] is a complex variable. For the 

fixed link 9 1 = 0 . 

Differentiating Equation 5.7 successively with respect to time, we get 

co 2 Z 2 + <0 3 Z 3 - co 4 Z 4 = 0 [5.8] 

and 

( a2 + i cq 2 2 )Z 2 + ( ct 3 + i ft)3 2 )Z 3 -( ou + i ®4 2 )Z4 = 0 
or 

V 2 Z 2 + V 3 z 3 - v 4 Z4 = 0 [5-9 ] 

where Vn [ = ( a n + i co n 2 ), n = 2, 3, 4 ] is a complex variable. 

From Equations 5.7 and 5.8, we get 


<y 2 4 Z 2 + <d 34 Z3 + 6) 4 Z| — 0 


[5.10] 
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and similarly, from Equations 5.7 and 5.9, we obtain 

V24 Z 2 + V34 Z3 + V 4 Zi = 0 [5-11] 

where co V] = 0 i x - co j and Vy = V; - Vj. 

Eliminating Z 2 , from Equations 5.10 and 5. 1 1, we have 

Z3 = [ { - V 4 + ( 6)4/ CO 24 ) V 24 } / { V34 - ( o) 34 / <y 24 ) V 24 } ] Zi. 

' [ 5 . 12 ] 

Then, from Equation 5.10 

Z 2 = - [ ( (0 4 / CO 24 ) Z( + ( Cl) 34 / (0 24 ) Z3 ] [5.13] 

and from Equation 5.7 

Z 4 = -Z,+Z 2 + Z 3 . [5.14] 

Thus knowing Z 2 , Z3 and Z4 ( as a function of Z\ ), we can draw the mechanism assum- 
ing a suitable value for one of the link-lengths. 

Flow chart with illustrations - The steps used in the programming are clear from Chart 
5.1. After computing the link-lengths, one can use the same software for its kinematic 
analysis for verifying the results. 

For illustration of case 1, we have taken the following data’s to design a mechanism. 

( 0 2 = 9O 0 , co 2 = 2 rad /sec and a 2 = 0rad/sec 2 ), 

( 04 = 60°, co 4 = 3 rad / sec and a 4 = -1 rad/ sec 2 ) 


and Li=10unit. 
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The results obtained, are 

L-2 = 26.602 unit, 

L3 = 23.394 unit 
and L4 = 13.654 unit ; 

which form a non Grashofian mechanism. The angular displacements and velocities can 
be verified from Plot 5.1 and Plot 5.2, respectively. 

For case 2, we have considered a problem with 

( co 2 = 1 0 rad / sec and a 2 = 40 rad / sec 2 ), 

( co 3 = 4 rad / sec and a 3 = 60 rad / sec 2 ), 

( co 4 = 20 rad / sec and a 2 = 0 rad / sec 2 ) 
and Li = 10unit. 

The results obtained, are 

L 2 = 18.256 unit, 

L 3 = 22.419 unit 
and L 4 = 5.2147 unit ; 
which give a mechanism of the crank-rocker type. 
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f PROBLEMS WITH KINEMATIC VARIABLES ] 



[Chart 5.1 ] 















[ Plot 5.1 ] 
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[ Plot 5.2 ] 



6. CONCLUSIONS AND 
SUGGESTIONS FOR FUTURE WORK 

6.1 CONCLUSIONS 

On the basis of the results presented, the following conclusions can be drawn. 

[ i ] The slider-crank mechanism for specified dead-center configurations under different 
eases ( discussed in section 2.1 ) can be synthesized and visualized with animation. 
Having all the three inputs ( H, Qrr and one of the three parameters E, L2 and L3 ), the 
software checks the existence of a mechanism. Also an optimum design for prescribed 
values H and Qrr can be obtained. 

[ ii ] The synthesis of a crank-rocker mechanism for dead-center problems has been done 
without using the Volmer’s nomogram. In comparison to the conventional nomogram 
approach, the new approach gives more accurate values of /3 op t, ( /i min )max and hence the 
link-lengths for better transmission quality. The accuracy can also be controlled at the 
cost of execution load of the program as number of iterations will increase. 

[ iii ] The overall optimization approach for the path generation problems can design 4R 
mechanisms for as many specified points as one wants, but the computational load will be 
increased. The present software also takes care of some constraints like the maximum- 
minimum link-lengths ratio and the compulsion for the designed mechanism to be a 

crank-rocker. 

In the present software, an option to weight the error function ( used for the least 
square technique ) according to the errors has been provided. In this way the error can be 
further reduced. If there is any critical point through which the coupler point must pass, 
this approach will be useful but the errors at the other points may be increased. 
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Sometimes it is possible that we can have more than one solution, in that case one 
can choose the suitable one iegarding other properties, like the transmission quality, the 
space 1 cquii oment, the sensitivity to the toleiance, the quick-return ratio etc.. If the 
designer gets a Ciiashotian double-rocket as a solution, then using the property of 
cognates one can have crank-rocker solutions for the same coupler curve. 

[ iv 1 Overall optimization approach for the function generation problem also gives good 
results, but it may require some trials on the initial values. Here, one can also visualize 
the nature of the function on the screen. If instead of a function, a set of coordinated 
inputs and outputs are given, then one can guess the type of the standard functions ( e.g. 
sinusoidal function ) by visualizing the curve ( 0 4 vs. 0 2 )• 

( v | Regarding the two cusps problem, one can definitely obtain a crank-rocker 
mechanism with the proper selection of X 02 and x B i ( section 4.1 ). We can also visualize 
the movement of the coupler point from one cusp to other on the screen itself. 

[ vi ] We can draw all the coupler curves of 4R and 3R-1P mechanisms. The coupler 
curves can be drawn in segments for some specified crank rotation, by which one can 
predict the relative velocity at different positions. 

[ vii ] The problems with kinematic constraints ( angular velocity, acceleration etc. ) can 
also be solved for the two cases discussed in chapter 5. 

6.2 SUGGESTIONS FOR FUTURE WORK 

[ i ] All the problems, which we have discussed for four-link mechanisms can arise in the 
four-link module of a six-link and other complex mechanisms. In that situation, by 
synthesizing the four-link module, one may be able to the synthesize the required 


mechanisms. 
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[ ii ] Besides the precision point and overall optimization approaches one can think about 
a third approach which includes the overall optimization and also takes care of some 
critical points. Some other constraints like the minimum transmission angle, tolerences 
etc. can also be included. 

[ iii ] The weighted error function can also be used in the case of function generation 
problem and again one can consider the constraints which are mentioned in the path 
generation problem. 

[ iv ] In the case of two cusps problem, since we can get more than one solution ( 
depending upon the assumptions X 02 and xbi ), like dead-center problems of crank-rocker 
mechanisms one can maximize the minimum transmission angle to get a better crank- 
rocker mechanism. 

[ v ] Like an atlas of coupler curves for crank-rocker mechanisms developed by Hrones 
and Nelson 1111 , one can think of storing all the coupler curves for 4R and 3R-1P 
mechanisms, which can be treated as a quick guide for the designers. 

[ vi ] One can design a 4R mechanism if kinematic conditions ( Angular displacement, 
velocity and acceleration ) for input and floating links are prescribed. This can be done by 
writing Freudenstein’s equation in terms of 02 and 03 . 



REFERENCES 


1. Mallik, A.K., Ghosh, A. and Dittrich, G., Kinematic Analysis and Synthesis of 
Mechanisms, CRC Press Inc., 1994. 

2. Deb, K., Optimization for Engineering Design, PHI, 1995. 

3. Norton, R. L., Design of Machinery, McGraw-Hill, Inc., 1992. 

4. Hartenberg, R.S. and Denavit, J., Kinematic Synthesis of linkages. New York, 1964. 

5. Hain, K., Applied Kinematics, McGraw-Hill, 1967. 

6. Sandor, G.N. and Erdman, A.G., Mechanism Design. Vol 1, Prentice Hall, 
Englewood Cliffs, New Jersey, 1991. 

7. Lai, J., Theory of Mechanisms and Machines, Metropolitan Book Co., 1992. 

8. Hall, A.S., Kinematics and Linkage Design, Prentice Hall, Englewood Cliffs, New 
Jersey, 1961. 

9. Tao, D.C., Applied Linkage Synthesis, Addison - Wesley Publishing Co., 1964. 

10. Lichtenheldt, W. and Luck, K., Konstruktionslehre der Getriebe ( in German ), 
Akademie-Verlag, 1979. 

11. Hrones, J.A. and Nelson, G.L., Analysis of the Four-Bar Linkages, MIT Press, 
Cambridge, Massachusetts, 1951. 



APPENDIX 


BASIC CLASSES OF C ++ USED IN THE SOFTWARE 

In developing a large and complicated software, classes are necessarily used. In' the 
present software, the basic classes used, are 

HI Revolute - It contains all the necessary characteristics of revolute pairs. 

The main member functions of this class are 

• EnterData ( ) ; This function asks for the link lengths through another class 
LinkLength. 

• IfExist ( ) ; It checks for the existence of a chain 

• FixLink ( ) ; After fixing a link chain becomes a mechanism. 

• MaxMin ( ) ; Finds the maximum and minimum lengths. 

• IsGrashof ( ) ; Checks for a grashofian mechanism. 

• CalAngle ( ) ; Calculates the output angle for each input angle, which is used 
in animation of 4R mechanism. 

• CalRange ( ) ; Calculates the range of the angular movement of any link of the 
mechanism. 

• EnterOmega ( ) ; Asks for the input angular velocity. 

• Enter Acc ( ) ; Asks for the input angular acceleration. 

There is an inherited class Prism of the class Revolute for prismatic pairs. 

[ii] LinkLength - It stores and returns the link length corresponding to a particular link. 

[iii] Range - It stores and returns the range of linear or angular movement of any link of 


the mechanism. 
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fiv] Animate - Using this class animation of the mechanisms along with creating of 
menus and windows are done. The initial interaction with the user is made through a 
member function of this class. The main functions of this class are 

• Draw_4R ( ) ; Animates 4R mechanisms. 

• Draw_3RlP ( ) ; Animates 3R-1P mechanisms. 

• Curve_4R ( ) ; Plots coupler curve for 4R mechanisms. 

• Curve_3RlP ( ) ; Plots coupler curve for 3R-1P mechanisms. 

• Cog4R ( ) ; Animates the cognates of crank-rocker and double-rocker 
mechanisms. 

• Display( ) ; Calculates the minimum transmission angle and display it along 
with the link lengths and types of mechanisms. 

• PlotTh4 ( ) ; It plots output angle vs input angles. 

• PIotVelGraph ( ) ; Plots input and output angular velocities corresponding to 
each input angular displacement. 

• GiveAcc ( ) ; It does acceleration analysis. 

[v] ScreenData - This class interacts with the user and asks for the data’s available 
through some member functions of class synthesis! to synthesis a mechanism for 
different types of problems, 

[vi] Angle - It stores and returns the angular data’s, which are used in synthesis of 
mechanisms. 

[vii] Point - It stores and returns the linear data’s or co-ordinates of points (required in 
path generation problems, cusp problems etc. ), which are used in synthesis of 
mechanisms. 

[viii] Synthesisl - This class is mainly used to compute the unknown design parameters 
lor all types of synthesis problems which we have already discussed in this work, the 
main functions of this class are 



|Mf # * 
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• AddPointData ( ) ; With the help of this, the required linear inputs for different 
types of problems are asked. 

• AddAngleData ( ) ; Same as above but this is for angular data ‘s. 

• Deadsol ( ) , Computes the unknown design parameters for dead centre 
problems. 

• SolveEq 1 ( ) ; For overall optimisation path generation problems 

• SolveEq2 ( ) ;. For overall optimisation function generation problems. 

• Cuspsol ( ) ; For two cusps problems. 

• . Kinsol ( ) ; For problems with kinematic constraints. 

Thus in the case of synthesis of mechanisms, ScreenData interacts with the user and 
executes the member functions of Synthesis] ( the member functions of Synthesis 1 
execute the functions of classes Angle and Point ) and transfer all the design parameters 
to class Animate for animation of the mechanism through EnterData ( ). 



